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Influence of corotation on the high energy synchrotron 
emission in Crab-like pulsars 
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ABSTRACT 

For Crab- like pulsars we consider the synchrotron mechanism influenced by relativistic 
effects of rotation to study the production of the very high energy (VHE) pulsed radi- 
ation. The process of quasi-linear diffusion (QLD) is applied to prevent the damping 
of the synchrotron emission due to extremely strong magnetic field. By examining the 
kinetic equation governing the QLD, apart from the synchrotron radiative force, we 
taken into account the the so-called reaction force, that is responsible for corotation 
and influences plasma processes in the nearby zone of the light cylinder (LC) surface. 
We have found that the relativistic effects of rotation significantly change efficiency of 
the quasi-linear diffusion. In particular, examining magnetospheric parameters typical 
for Crab-like pulsars, it has been shown that unlike the situation, where relativis- 
tic effects of rotation are not important, on the LC surface, the relativistic electrons 
via the synchrotron mechanism may produce photons even in the TeV domain. It is 
shown that the VHE radiation is strongly correlated with the relatively low frequency 
emission. 

Key words: Pulsars: general - Acceleration of particles - Radiation mechanisms: 
non-thermal. 
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1 INTRODUCTION 

According to the standard pulsar model, in particular, the 
so-called polar cap model |Sturrock 197ll ). particles uproot 
from the st ar's surface, and accelerate alon g the magnetic 
field lines jRuderman fc Sutherland TOTS ') inside a gap 
zone. In general, it is believed that these particles produce 
the nonthermal radiation, w hich is interpreted in terms o f 
the synchrotron mechanism JPacini 197ll: IShklovskv 197C ) 
and the Inverse Compton Scattering (jBlandford et al. 199QT ) 
respectively. Generally speaking, due to strong synchrotron 
losses, relativistic electrons quickly lose their perpendicular 
momentum, and transit to their ground Landau states. 
Therefore, particles may be approximately described as 
moving one-dimensionally along the field lines, which in 
turn means that the synchrotron mechanism does not 
contribute in high energy radiation processes. On the other 
hand, under certain conditions, due to the quasi-linear 
diffusion (QLD), the pitch angles might be recreated, lead- 
ing to the efficient synchrotron e mission. This method has 



already been applied to pul sars ( Machabe li fc Usov 197S 
Malov fc Machabeli 200ll; IChkheidze fc Machabeli 2007 
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and active ga lactic nuclei (|Osmanov fc M achabel i 2010l : 
lOsmanov 2010 ) . 

In the context of QLD, the r ecent detection o f the VHE 
radiation from the Crab pulsar (|Aliu et al. 20081 ) could be 
very important. MAGIC Cherenkov telescope from 2007 Oc- 
tober to 2008 February, has discovered the pulsed emis- 
sion above 25GeV, which reveals several characteristic fea- 
tures. A special interest deserves the coincidence of signals 
from different ener gy bands rangin g from radio to VHE 
(> 25GeV) domain (|Aliu et al. 20081 ). The authors conclude 
that according to the data, the polar cap models must be 
excluded from the possible scenario of the radiation, which 
must happen far out in the magnetosphere. Analysis of the 
MAGIC data also implies that a location of the aforemen- 
tioned VHE and low energy emissions must be the same. 
By taking into account the synchrotron reaction force, we 
constructed the kinetic equation governing the QLD and the 
results of the MAGIC data have been interpreted (Macha- 
beli & Osmanov 2009). We argued that the observed VHE 
radiation is produced by the synchrotron mechanism, having 
properties, which are in a good agreement with the obser- 
vations. The following work (Machabeli & Osmanov 2010) 
was related to the problem of the curvature radiation and 
the inverse Compton scattering in the context of the same 
observations. We have shown that these mechanisms do not 
contribute to the VHE domain detected by MAGIC. 
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In general, magnetic field in pulsar magneto-spheres is 
huge and ranges from I0 6 G (far out in the magnetosphere) 
to 10 12 G close to the neutron star's surface. Therefore, mag- 
netospheric plasma is in the frozen-in condition and is forced 
to follow the magnetic field lines. This means that effects 
of corotation in plasma acceleration and emission processes 
could be very important. In general, it is believed that the 
high energy radiation comes from an area located relatively 
close to the LC surface, where corotation is extremely sig- 
nificant. For this purpose it is reasonable to consider the 
quasi-linear diffusion by taking relativistic effects of rota- 
tion into account and see how the corotation influences the 
mentioned processes. 

In the present paper we study the role of the corota- 
tion in the quasi-linear diffusion and we show that in the 
very vicinity of the LC surface, contrary to the quasi-linear 
diffusion, effects of corotation attempt to decrease the pitch 
angles, and kill the subsequent synchrotron emission. 

The paper is organized as follows. In Section 2 we con- 
sider the quasi-linear diffusion by taking the effect of coro- 
tation into account, in Sect. 3 we present our results and in 
Sect. 4 we summarize them. 



2 MAIN CONSIDERATION 

In this section we present the model, generalizing our previ- 
ous approach by taking the fact of corotation in to account 
and see how the efficiency of QLD depends on it. 

Generally speaking, the pulsar magnetospheres are 
composed of low and high energy particles respectively. 
Therefore, in the framework of the paper for simplicity we 
assume that the magnetosphere is consist of two compo- 
nents: (a) the plasma component with the Lorentz factor, 
7 P and (b) the beam component with the Lorentz factor, i b . 

According to the mechanism of the quasi-linear diffusion 
the following transverse mode generates 



ut w kc (1 — S) . 



8 = 



(1) 



where k denotes the modulus of the wave vector, c is the 
speed of light, lo p = Ann p e 2 /m is the plasma frequency, e 
and m are electron's charge and the rest mass, respectively, 
n v is the plasma density, lob = eB/mc is the cyclotron fre- 
quency and B is the magnetic induction. It can be shown 
that the aforementioned wa ves excite if the cyc lotron reso- 
nance condition is satisfied (|Kazbegi et al. 1992) 



fen v - 



k x u ± — — = 0. 

lb 



(2) 



By and k± we denote the wave vector's longitudinal (par- 
allel to the background magnetic field) and transverse (per- 
pendicular to the background magnetic field) components 
respectively, u = cvib/ Plob is the drift velocity of elec- 
trons, v is the component of velocity along the magnetic 
field lines and p is field line's curvature radius. In this pa- 
per we consider the beam particles as to be the resonance 
particles. One can show from Eqs. (1 1 121) that the excited 
cyclotron wave is characte rized by the following frequency 
l|Malov fc Machabeli 200lh 
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(3) 



Generally speaking, particles moving in the magnetic 
field, experience two forces, one of which, G , having the 
following components ()Landau &: Lifshitz 197ll ) 
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is responsible for the conservation of the adiabatic invari- 
ant, I = 3cp 2 / 2eB and the second one is the synchrotron 
radiative force l|Landau fc Lifshitz 197ll ) 
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where a = 2e 2 tjg/(3c 2 ). But one can show that for the 
Crab- like pulsars the radiation reaction force, |F|, exceeds 
|G| by many orders of magnitude (Machabeli & Osmanov 
2009,2010). 

On the other hand, apart from the synchrotron radiative 
force, relativistic electrons also experience a reaction force 
responsible for the corotatiorQ This force can be estimated 
by a simple mechanical analogy intr oduced by Machabe li & 
Rogava (1994) and reconsidered by (|Rogava et al. 20031 ). In 
the framework of this approach, instead of magnetic field 
lines and particles sliding along them, one considers a coro- 
tating pipe with a bead inside it. If we assume that the 
magnetic field lines are straight, then the reaction force act- 
ing on a single particle from the beam component can be 
given by (|Rogava et al. 20031 ) 



R 



dp ± 
dt 



(6) 



where Q is the angular velocity of rotation, r is the radial 
coordinate, 



p n = ~/ b mv 
and 

p ± = ibmrQ 



(7) 



(8) 



are momentum's longitudinal and transversal components 
respectively. 

Considering a single particle approach, one can show 
that the Lorentz factor of a particle moving along the 
corotating straight magn etic field yields the equation 



ijRieeer fc Mannheim 200ol ) 



lb = 



where 



(9) 



c 



ro and vq are particle's initial position and the initial ra- 
dial velocity respectively and r; c is the light cylinder radius. 
From Eq. ([9| we see that the closer to the LC, the bigger 
the Lorentz factor. Therefore a dynamical effect of the coro- 
tation becomes extremely efficient nearby this area [see Eq. 

] where the reaction force can be approximated as to be 
[sec Eq. ()A6|I in Appendix] 



1 In the local frame of reference particles experience the centrifu- 
gal force, whereas considering dynamics in the laboratory frame, 
one has to examine the reaction force acting on the electrons from 
the magnetic field line. 
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Figure 1. The synchrotron emission energy versus the Lorcntz 
factors of beam electrons. The set of parameters is P ss 0.0332s, 
B ss 1.7 x 10 7 G, 7t,o = {1, 10 2 , 10 4 }. 
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As we have mentioned, in the Crab pulsar magneto- 
sphere one may neglect the effects of G in comparison with 
F. Therefore, the process of the quasi- linear diffusion is 
mainly influenced by the synchrotron radiative force, F and 
the reaction force, R. The synchrotron radiation reaction 
force, acting on a particle, attempts to decelerate it, there- 
fore relativistic electrons loose their perpendicular momen- 
tum, leading to an inevitable decrease of the pitch angles. 
In spite of the fact that contrary to F, the reaction force 
plays an accelerative role in electron's dynamics, in the con- 
text of the pitch angles its role is the same. In particular, 
as we have already discussed, the mentioned reaction force 
is a direct consequence of corotation. This means that R is 
constructed in such a way that particles always must stay 
on the magnetic field lines, therefore, the reaction force also 
attempts to decrease the pitch angles. On the other hand, 
the quasi-linear diffusion, that arises through the influence 
of the waves back on the particles, via the QLD attempts 
to widen a range of the pitch angles. The dynamical process 
saturates when the effects of the above mentioned forces 
are balanced by the diffusion. For 7^ 3> 1 it is easy to 
show that in the quasi-stationary case (d/dt = 0), the cor- 
respond ing kinetic equation gove rning the QLD writes as 
follows l|Malov fc Machabeli 2o"oih 



1 d[iP(R + F ± ) X ] = d_ 
p.. dip dip 



(11) 



where x = x(V0 i s the distribution function of particles with 
respect to the pitch angles, 



7re ni,c 



(12) 



is the diffusion coefficient, nt = B/ Pee is the density of the 
beam component and P is the pulsar's rotation period. The 
solution of Eq. (1121) writes as follows 



mc~f b R 



As is clear from this expression, electrons are differently 
distributed for different pitch angles. Let us note that the 
reaction force and the synchrotron radiative force have op- 
posite directions, since Fj_ is responsible for deceleration of 
electrons, whereas R accelerates them. It is clear that close 
to the LC, a value of \R\ oc j^ 2 is very high and for very 
small pitch angles exceeds \F±\ and therefore, higher val- 
ues of pitch angles lead to higher values of the distribution 
function. But by increasing ip the corresponding reaction 
force does not change [see Eq. ([10}], whereas the synchrotron 
radiative force is very sensitive to this change, F± ~ ip 3 . 
Therefore, at a certain pitch angle, ipo, these forces will bal- 
ance each other and the distribution function will reach its 
maximum value. In particular, one can straightforwardly 
show that the expression of the distribution function [see 
Eq. (JT3J] peaks at 



1 



2mcSl 

2 V2 
"76076 



1/3 



(14) 



which indeed corresponds to R — F±. Here, we have taken 
into account the following relation m sa l/7 2 [ro/ri c <S 1, 
see Eq. ©]. 

It is reasonable to estimate how steep is the distribu- 
tion function. For simplicity we consider the case 760 = 1, 
then examining the Crab pulsar's (P sa 0.0332s) magneto- 
spheric parameters close to the LC surface: B S3 1.7 x 10 7 G, 
76 = 10 7 , one can see that at the peak, ipo S3 0.02rad, 
the distribution function exceeds its value at ip = by 
many orders of magnitude. This means that most of the 
particles will be characterized by the peak value of the 
pitch angle, and therefore, its average value can be es- 
timated as to be ipo. We have taken into account that 
B = B a tRl t /r't, where R Bt « 10 6 cm is the pulsar's radius, 

B st S3 1.8 x 10 12 v P x P15G is the magnetic induction on 
the star's s urface and As = 10 15 P (A s = 421 for the Crab 
pulsar, see ^Manchester fc Taylor 198Ch ). 

Since we consider centrifugally accelerated particles, our 
approach is based on the assumption that electrons coro- 
tate with the spinning pulsar, which is valid only inside a 
certain zone (corotation zone with a radius r c ), because in 
this region the magnetic field is strong enough to channel 
the flow. This means that, our approach is valid if the fol- 
lowing condition is satisfied: (r; c — r c )/r; c <C 1. It is clear 
that corotation takes place if the magnetic energy density 
exceeds the beam energy density, therefore, r c can be can 
be estimated by the condition: B(r) 2 /(8n) ss mc 2 nf > y(r). If 
we consider a range of 760 = (1; 10 4 ), then by taking Eq. 
((9)1 into account, one can see that for 76 ~ 10 r one has 
(rjc — T c )/ri c ~ (10 -11 ; 10 -7 ) <C 1, which means that coro- 
tation is violated in a very close region to the LC surface and 
for almost the whole course of motion the spinning magnetic 
field channels the electrons. 

When relativistic particles move in the mag- 
netic field, they emit electromagnetic waves, which in 
our case correspond to the following photon energies 
l|Rvbicki fc Lightman 1979h : 



e Te v ~ 1-7 x 10^ 20 B7 2 sin^ 



v RV3 v „,ll/6 x 7 -2/3_ 



10- 16 x B 1/3 x 7 11 



(15) 
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Figure 2. The synchrotron emission energy versus the cyclotron 
frequency. The set of parameters is P ss 0.0332s, B rs 1.7 X 10 7 G, 
760 = {1,10 2 ,10 4 }. 



3 DISCUSSION 

In this section we are going to apply the method to the 
Crab- like pulsars and see what changes in the process of the 
QLD when corotation is taken into account. 

In Fig. [T] we show the dependence of photon energies on 
the Lorentz factors of the beam electrons for three different 
initial values of the beam Lorentz factors. The set of param- 
eters is P « 0.0332s, B « 1.7 X 10 7 G, -y b0 = {1,10 2 ,10 4 }. 
As is clear from the plots by increasing 7(, the corresponding 
emission energy increases as well. This is a natural result, 
because more energetic electrons emit more energetic pho- 
tons [see Eq. (|15[) ]. A second feature concerns the behaviour 
of erev versus 760- It is evident that the higher the initial 
beam Lorentz factor the lower the emission energy. Indeed, 
as one can see from Eq. (|14[) the pitch angle behaves as 
7~ 2 / 3 , therefore, for higher initial beam Lorentz factors one 
obtains lower values of the pitch angles and hence, the corre- 
sponding synchrotron energy becomes lower. The relativistic 
effects of rotation could also be interesting in the context of 
the recently detected VHE pulsed emission (> 25GeV) from 
the Crab pulsar. According to the interpretation presented 
in (Machabeli & Osmanov 2009,2010), the mentioned radi- 
ation is formed on lengthscales of the order of 10 s cm via 
the quasi-linear diffusion. As it has been shown, in order 
to explain the data, the value of the beam Lorentz factor 
must be of the order of 3 x 10 8 . For centrifugally accelerated 
electrons the effects of corotation become significant in the 
nearby zone of the LC surface. In particular, from Eq. (|10l) it 
is clear that the reaction force is proportional to "f^ 2 , which 
in turn asymptotically increases by reaching the LC surface 
[see Eq. The corresponding lengthscale of a layer on 

the light cylinder zone, where the effects are intensified by 
corotation, can be estimated as follows 



A 



76 



I'lc 



dj b /dr " 2 7i) ' (16) 

where we have taken into account Eq. ((9}. As we see from 
this expression, the layer is quite thin, A <C r; c . Unlike the 
papers by Machabeli & Osmanov (2009,2010) where the ef- 



fects of corotation are not significant, and the VHE radiation 
(> 25GeV) can be achieved by the relativistic particles with 
76 ~ 3 x 10 s , a part of the magnetosphere located very close 
to the LC surface, may guarantee the aforementioned emis- 
sion energy for lower Lorentz factors. In particular, from 
Fig. Q] it is clear that the synchrotron process supported 
by the QLD produces 25GeV photons for 76 ~ 5.3 x 10 6 
(7i.o = 1), 7fc ~ 8 x 10 6 ( 7 6o = 10 2 ) and 7b ~ 1.2 x 10 7 
(760 = 10 4 ), which are much lower than 3 x 10 8 . The quasi- 
linear diffusion is so efficient that under certain conditions 
it can provide TeV photons as well. As one can see from the 
figure, the relativistic electrons with Lorentz factors 5 x 10 7 
(760 = 1) may provide TeV emission (1.5TeV). 

According to the present method, the process of the 
QLD depends on excitation of the unstable cyclotron waves, 
which correlate with the high energy radiation. In Fig. [2] we 
show the behaviour of synchrotron emission energy versus 
the cyclotron frequency. The set of parameters is the same 
as in Fig. [T] As is seen from the plots, the VHE radiation 
is strongly connected with relatively low energy radiation 
starting from microwave (5 x 10 10 Hz) to optical (10 14 Hz) 
domains respectively. 

The investigation shows that very close to the LC area, 
where effects of corotation are extremely important, the 
quasi-linear diffusion becomes so efficient that under favor- 
able conditions it may provide VHE radiation for relatively 
lower energy electrons than in the regime considered by 
Machabeli & Osmanov (2009,2010). The aim of the present 
paper was to demonstrate for Crab-like pulsars that the rel- 
ativistic effects of rotation are of fundamental importance 
for studying the VHE emission via QLD. Studying rotation- 
ally driven quasi-linear diffusion for 1-sec pulsars also is an 
interesting problem, but it is out of the scope of the present 
paper and sooner or later we are going to consider it as well. 



4 SUMMARY 

(i) We have examined the VHE radiation of Crab-like pul- 
sars via the quasi-linear diffusion intensified by the effects of 
corotation which takes place in the nearby zone of the LC 
surface. 

(ii) It has been emphasized that due to very strong mag- 
netic field of pulsar magnetospheres the efficient energy 
losses lead to the damping of the synchrotron process. Close 
to the LC surface the effects of corotation become signifi- 
cant and apart from the synchrotron radiative force, which 
decreases the pitch angles, also the reaction force (responsi- 
ble for corotation) has to be taken into account. Generaliz- 
ing the kinetic equation governing the QLD, we have found 
the particle distribution with respect to ip and estimated its 
average value. 

(iii) Considering Crab-like pulsar's magnetospheric pa- 
rameters we have found that the quasi-linear diffusion be- 
comes more efficient on the LC zone, than for locations far 
from this area. It has been shown that the synchrotron emis- 
sion process, via the QLD, might provide VHE emission even 
in the TeV domain. 

(iv) We have found that the VHE radiation is strongly 
connected with the cyclotron emission having relatively low 
frequencies. 
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APPENDIX A: 

In this section we would like to estimate the reaction force 
close to the LC surface. 

From Eqs. (08]) one can see that 

R = mi] ^] 11 ^ + mj b vQ, = mflr^- + 2my b vtt. (Al) 

By applying Eq. ([9]) for the nearby zone of the LC area, 
r ~ ric, the aforementioned expression writes as follows 

R « 2mrh 1,2 vQ, 1 l + 2m'y b vQ,, (A2) 

which by taking into account a natural relation m 1/,2 7& 3> 1 
reduces to 



R m 2mm 1 vfl-yi . 



(A3) 



Machabeli & Rogava (1994) showed (see Eq.10 in the corre- 
sponding article) 




(A4) 



which, applying to the LC surface, by using Eq. ((9]), obtains 
the following approximate form 



(A5) 




~1/4 1/2' 

m 1 / 4 7 () 

and reduces Eq. ()A3[) to the final expression 
R ~ 2mm 1 ^ A cQ.^ 2 . 
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